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Localization is a phenomenon long studied in condensed-matter physics [[T]]. Concerned with 
conduction in disordered media, it is based on the study of eigenstates of the Schrodinger equation 
in a random potential. Its counterpart in lattice gauge theory is the study of the spectrum and 
eigenstates of a Hermitian fermion kernel in the fluctuating gauge field. Here we deal with Wilson 
fermions, with kernel Hw = ysDw, in an assortment of gauge ensembles We find low-lying 

localized states, extending up to an energy called the mobility edge; and extended states at higher 
energy. We characterize the localized states by their localization length and their support length 
(related to their inverse participation ratio, or IPR). 

The Wilson fermion kernel is not so popular these days for its own sake but rather for its role in 
constructing domain- wall fermions [Q] and the closely related overlap fermions [|5]]. The simplest 
overlap kernel is given by 

^ov = l-75Sgn(//w)- (1) 
The range of Dqv can become long if the spectrum of Hw contains low-lying (or zero) modes |^]. 
We will argue here that, though the energy threshhold of Hy^ lie at zero, the range of Dqv is deter- 
mined by the position of the mobility edge — that is, the low-lying modes are generally harmless. 
The mobility edge perforce lies at non-zero energy, as long as one avoids the Aoki phase [^]. In- 
deed, the descent of the mobility edge to zero serves as a useful definition of the onset of the Aoki 
phase. 

1. Localization basics 

As stated above, we address the problem of the spectrum of H^/ in the random gauge field 
Ux^- Dislocations in any given configuration Uxp. can create bound states, even zero modes [||], in 
the spectrum of Hv/. Of course we do not study single configurations of Uxfj. but averages over an 
ensemble. The result is localization. 

For a general overview, let's look at Fig. |l[ We assume that the free operator has a gap for 
the parameters we choose; above that gap lie plane-wave states. Condensed-matter physicists call 
the top of the gap the band edge. If we add a single dislocation to the otherwise constant gauge 
field, the plane waves will be replaced by scattering states and a bound state might appear in the 
gap. A large number of dislocations will create a large number of bound states, and the scattering 
states will show effects of multiple scattering. In finite volume, the scattering states as well as the 
bound states form a discrete spectrum. 

In the infinite volume limit, with a fixed density of dislocations, both the bound states and the 
scattering states form a continuum. The energy that separates the two is called the mobility edge. 
The infinite volume limit automatically gives us an average over the shapes and positions of the 
dislocations in any finite subvolume. Alternatively, we can keep the volume finite and average over 
gauge field configurations ourselves, as is done in lattice gauge theory. 

H\f can only have zero eigenvalues in the supercritical region, — 8 < mo < 0. As it happens, 
this is just the region of interest for the construction of domain-wall and overlap fermions. The 
Aoki phase lies in this supercritical region; we stay outside the Aoki phase so that the mobility 
edge is above zero. We will present numbers below for mo = —1.5. For the gauge couplings 
we choose, the theory lies between the Aoki "fingers." These couplings are in fact popular among 
users of domain-wall and overlap fermions. 
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Figure 1: How the spectrum of Hw changes as one adds (top to bottom) one dislocation, many dislocations, 
and a random ensemble of dislocations to the gauge field. 



2. The spectral density 

Since averaging over the gauge field smears the energy eigenvalues, there is never a discrete 
spectrum of bound states accompanied by a continuum of extended states. The spectrum in the 
disordered system is described by a continuous density of states, 

pW = Ul,5{X-Xn)), (2.1) 

n 

where A is the energy eigenvalue. We show in Fig. § the density of states for the gauge ensemble 
generated with the Wilson plaquette action. The vertical bars indicate the mobility edge for the two 
couplings. Nothing special happens at the mobility edge. 
[We find it convenient to calculate p via the identity 

7rp(A) = ^ lim (ImTrG(A + (2.2) 

V HI 1^0 

where G{z) = {Hw — z)^ is the resolvent of Hw- The mobility edge is where the averaged local- 
ization length /^(A) (see below) diverges.] 

The main point of Fig. ^ is to show that there is indeed a nonzero density of states at zero 
energy. Table |l| shows how this density changes with cutoff and with the choice of gauge action. 
Improved actions dramatically lower p(0), but in no case is it actually zero. 

In a theory with dynamical Wilson fermions, these zero modes will be of no consequence 
because they will cause the fermion determinant to vanish. In any other theory, however, the 
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Figure 2: Eigenvalue density p{X), in lattice units, for j3 = 5.7 (a ' = 1 GeV, upper curve) and j3 = 6.0 
(fl^^ = 2 GeV, lower curve). Plaquette action. 



Table 1: Eigenvalue density at A = for three gauge actions, each at two values of the gauge coupling. 

spectrum of H]/v does not directly determine the fermion determinant. Whether the ensemble is 
quenched, or whether it contains the domain-wall or overlap fermion determinant, the ensemble 
average will furnish H]/v with modes near (and at) zero energy. Will this destroy the locality of 
Dov? No — because (outside the Aoki phase) these modes are localized. 

3. To the overlap kernel 

A localized mode can be described by its support length Ig, which contains most of the mode's 
density, and by its localization length Ig, which is the decay length of any long-range tail (see 
Fig. 111). Each can be averaged over all eigenmodes with eigenvalue A, giving k{X) and /^(A). The 
average support length can be used to determine in what range of A the localized modes are dilute. 
In view of the mode density shown in Fig. ^, it is clear that low-lying modes will be dilute while 
the higher modes will be more crowded. 

The dilute, the less dilute, and the extended modes all contribute to the overlap kernel Hov 
The contribution of the dilute modes, those below some energy A, is of the form 



Wilson Iwasaki DBW2 



a 



a 



IGeV: p(0) = 10-2 ^ 6-10-^ IQ-^ 
2 GeV: p(0) = 10-4 7 • 10"^ < 10 





(3.1) 
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Figure 3: Characteristic lengths of a locaHzed eigenfunction: the support length l^, which is the size of the 
region containing most of the amplitude; and the localization length l(, which is the decay length of the tail. 

where the second line follows from the steep rise in p(A) seen in Fig. ^. Note that these isolated 
modes contribute an exponential tail characterized by their localization length; their low energy is 
harmless.^ The dense modes are up near the mobility edge Xc, if we lump them with the continuum 
of extended modes then we can estimate that the fall-off of their contribution will be governed by 
Xc, viz. 

(|Dov(jc,3')|)ext ~'^exp(-Af|x-j|), (3.2) 

where ^ is 0(1). 

Let's put in numbers for the plaquette action at j8 = 6.0. We measure the mobility edge to be 
at Ac 0.41 and we fix the demarcation point to be A 0.2. The result is 

(|Dov(x,3')|) ^ 10-^exp (-^) + ^(1) -exp (-^) • (3-3) 

The second term, governed by the mobility edge, wins — both in prefactor and in exponent. We find 
this to be true in all the cases (actions and couplings) we have studied, including those in Table 
1. The reverse could be true in other cases; then the range of Dqv would be determined by the 
representative localization length, le(X). The fact that the gap is zero would still be harmless. 

We can interpret Xc as a mass scale, the mass of effective excitations that influence the overlap 
kernel. When the cutoff is =2 GeV, we find that Xca^^ ~ 800 MeV for all three actions. It 
may be a disappointment to find that the cutoff is felt at an energy that is so low, but this is still well 
above Aqcd which is the true dynamical scale. One ought to worry, however, at stronger couplings. 
When a^^ = I GeV, the value of Xc-a^^ comes out to be only 250-320 MeV, not a place we would 
like to see unphysical particles in the spectrum. 

4. Conclusions 

The three main lessons of our work: 

' The contribution of dilute localized modes to any Green function will decay in space with the localization length 
Ig. This is in fact the most convenient way to measure by calculating almost any Green function P]. 
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1. The mobility edge at Ac > [or /^(A), if it is greater than ^] assures a finite range for Z)ov> 
even though has no gap in the disordered gauge field. 

2. One should demand that Xca~^ » -'^QCD- 

3. XcOT^ is fairly insensitive to the gauge action, even as p(0) varies widely. 

For discussions of these issues in the context of current large-scale lattice computations, see the 
contributions of T. Draper and of P. Boyle to these proceedings. 
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